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^ ! Abstract 

Recently there has been a proposal for modified gravitational /(it!) actions which 
include a direct coupling between the matter action and the Ricci scalar, R. Of 
particular interest is the specific case where both the action and the coupling are 
Q^i linear in R. It is shown that such an action leads to a theory of gravity which includes 

^ ' higher order derivatives of the matter fields without introducing more dynamics in 

, the gravity sector and, therefore, cannot be a viable theory for gravitation. 
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Q ' Even thought alternative theories of gravity have always been a subject of 

. study since the introduction of General Relativity (GR), nowadays they ap- 

pear to attract increased interest. This is mainly due to the fact that they 
appear as candidates for explaining the nature of dark matter and dark en- 
^ ■ ergy and, consequently, they might bring us closer to answering some of the 

c3 ' most important open questions of late time cosmology. 

One of the proposals which have been given much attention recently is that 
of /(i?) theories of gravity: theories described by the action 

Sf = I d^xV^[^f{R) + L„,{g,,, ^)] (1) 

where g is the determinant of the metric g^^, k = 8ttG, f{R) is a general 
function of the the Ricci scalar R, is the matter action and ip collectively 
denotes the matter fields. There exist more than one variational principle 
which could be applied to the above action. Namely, one could use the standard 
metric variation, in which case the result is the so-called metric f{R) theories 
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of gravity, or an independent variation of the metric and the affine connection 
(Palatini variation) which leads to Palatini f{R) gravity (see [1,2,3,4,5] for 
early works and review . In the Palatini formalism the connection F^^^, 
is independent from the metric. Therefore, R is replaced in action (1) by 
TZ = g^'^lZ^y, where TZfj^u is the Ricci tensor constructed with the independent 
connection. Note that, in action (1), Lm is taken a priori to be independent 
of the connection. If this assumption is abandoned the outcome is yet another 
version of f{R) theory, metric-affine f{R) gravity [7]. Needless to say that all 
three version of f{R) gravity are different theories, despite the similarity of 
their actions [8,9]. 

Much work has been done in studying the phenomenology and the viability 
of f{R) gravity (see [10] for some examples). Even though it is not very likely 
that some action within this class can constitute a complete and viable alter- 
native to GR, f{R) theories have proved to be very useful as toy-theories of 
gravity. As such, they can significantly contribute to out understanding of the 
gravitational interaction. 

Recently, the f{R) idea was generalized to include a coupling between the 
Ricci scalar and the matter itself. The action presented in [11] was 



where, /i^2 are arbitrary functions of the R, L.^ is the matter action and A 
acts as a coupling constant. Other actions with similar couplings between R 
and matter have been studied in [12,13,14]. 

Variation of the action (2) with respect to the metric yields 



is the matter stress-energy tensor as usual. Note that T^jy is not divergence- 
free here, due to the direct coupling between matter and R. Therefore, the 
theory will exhibit violations of the Einstein Equivalence Principle [15], which, 
however, can in principle be controlled by tuning A [11]. 

^ It is worth pointing out that f{R) gravity in the metric formalism appears to be 
unique among higher order gravity theories in that it avoids the so-called Ostro- 
gradski instabiUty [6]. 




(2) 



F,{R)R^, - -h{R)g^u - V^V,Fi(i?) + g^^uF^iR) 

= -2XF2{R)L^R^, + 2A(V^V^ - g^,U)L^F2{R) 

+ [l + A/2(i?)]T^., 



(3) 
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In [11] the authors examined whether theories described by the action (2) 
could account for phenomenology usually attributed to dark matter. In [16] 
the consequences of such a theory in stellar equilibrium were studied. As a 
matter of fact, the action used there was a particular case of action (2): the 
one in whicn ^ I fi{R) = k~^R , and f2{R) = R- It was argued in [16] that this 
linear choice for /i and /2 should serve as a good approximation, since, firstly 
it was expected that the phenomenology related to the coupling should over- 
whelm any phenomenology coming from the non-linearity of /i and, secondly. 
Solar System and other equivalence principle related tests will force A/2 -C 1, 
which combined with the first requirement yields f2{R) = -R as the simplest 
possibility. 

With these choices for /i 2 the action reads 

S = j ^gd^x 1^ + [1 + Ai?] L„| , (5) 
and the corresponding field equations are 

+ [1 + Ai?]/€T^„ (6) 

Let us first perform a straightforward comparison between eq. (3) and eq. (6). 
When fi is linear Fi is a constant and, therefore, the last two terms in the left 
hand side (Ihs) of eq. (3) are no longer present in eq. (6). Additionally, when 
/2 is linear, F2 is constant as well, so R is no longer present in the second term 
on the right hand side (rhs) of eq. (6). 

What is important to notice, however, is that these are the terms which con- 
tained derivatives of R. Indeed, eq. (3) is fourth-order in the metric, containing 
second derivatives of R. On the contrary, eq. (6) is just second order in the 
metric, exactly like GR. 

Eq. (3) differs from Einstein's equation in another important way as well: it 
contains higher derivative of the matter fields through the second term in the 
rhs. This characteristic still remain even when /i 2 are chosen to be linear, as 
can easily be seen in eq. (6). Assuming, as usual, that the matter action Lm 
contains only to first derivatives of the matter fields (so as to lead to second 
order field equations for the matter fields when varied with respect to them), 
eq. (3) is a fourth order partial differential equation (PDE) in the metric and 
a third order PDE in the matter fields. Eq. (6) is again a third order PDE in 
the matter fields but just a second order PDE in the metric. As a reference 



^ We are using here the conventional notation for k, whereas in [11,16] k is taken 
to be (167rG)-i. 
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point, Einstein equations are second order PDEs in the metric and first order 
in the matter fields. 



A serious issue becomes immediately apparent for the case where /i 2 are 
linear: the differential order in the matter fields is higher than the differen- 
tial order in the metric! Remarkably this has gone unnoticed until now, even 
though it is a source of serious problems. 

To make this clearer, let us recall some past experience form Palatini f{R) 
gravity. Applying the Palatini variational principle on action (1) one gets 

f'{n)R^,,) - ^f{n)g^, = «:T^„ (7) 



^f'{n)g^'') = 0, (8) 

where Va is the covariant derivative defined with the independent connection. 
It is easy to show (see for instance [5,17]) that these equations can actually 
be combined to a single one: 

- Ijr2 ((V,/')(V./') - Ig.A'^fr) , (9) 

where, however, 71 and consequently f(jV) and f'(jV) are actually functions 
of T = g^'^T^^, due to the fact that the trace of eq. (7) 

/'(7^)i^-2/(7^) = /^T, (10) 

relates TZ and T algebraically (for instance, for f{TZ) = TZ + eTZ"^, where e is 
some parameter with suitable dimensions, R = —kT). Therefore, the Ihs of 
eq. (9) depends only on the matter fields and at the same time includes, up 
to third derivatives, exactly like the Ihs of eq. (6). At the same time both of 
these equation are second order PDEs in the metric. 

Now, the presence of the higher order derivatives of matter fields in eq. (9) 
has been shown to lead to very serious shortcomings for Palatini f{R) gravity. 
The post-Newtonian metric becomes algebraically related to the density [18] 
making the result of Solar System tests density dependent. This is clearly un- 
acceptable since such test should be valid over a wide range of densitieiil. At 
the same time, the presence of these higher order derivatives has been shown 
to lead to the appearance of singularities on the surface of specific spherically 
symmetric matter configurations when one attempts to match an interior to 
the unique exterior solution, making any matching impossible [17]. This leads 



^ In [19] concerns have been expressed on whether issues of averaging can have an 
affect on this but it is still questionable if this is the case. 
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to the absence of solution describing the gravitational field of physical objects 
which can even be described with Newtonian gravity [17]. Finally, even the 
problem found earlier in Palatini f{R), the unacceptable modification it ap- 
pears to introduce in the standard model of particle physics [20], seems to 
have its root in the presence of the higher order derivatives of matter fields 
[17]. 

It has to be stressed that these problems are not specific to Palatini f{R) grav- 
ity (or to the dynamically equivalent ujq — —3/2 Brans-Dicke theory [8]), but 
as already mentioned in [17], they are generic problem for theories in which 
the differential order of the matter fields in the field equations is higher than 
the differential order in the metric. This is because, the unusual behaviour 
discussed here does not really have to do with the fine details of the field 
equations, but merely with the fact that this unorthodox differential structure 
makes gravity non-commulative: the metric ceases to be an integral over the 
source and becomes algebraically related to the matter fields and their deriva- 
tives. Therefore, any discontinuities in the latter, which are in general allowed, 
will become unacceptable discontinuities and singularities in the metric, lead- 
ing to unacceptable phenomenology. 

Having said these, it is straightforward to realize that the theory described 
by the action (5) and the field equations (6) is bound to be burdened with 
essentially the same problems and no further calculation is actually needed to 
show that. Just for the sake of clarity and in order to further verify our claim, 
let us consider the weak gravity regime. 

First of all, notice that if one takes the trace of eq. (6) this gives 

_ GnXnLm - kT 

1 + kX{T-2LJ- ^''^ 

Replacing this back in eq. (6) in order to eliminate R yields after some ma- 
nipulations 



R 



1 + 2kAL^ 



+ 1 + kX{T - Lm) \ ~ 2^'" 



(12) 



This is, in general, the most usual and convenient form to discussed the weak 
field regime. Now, suppose the metric is taken to be flat plus a perturbation: 

Qnv = Vu-v + (13) 

We will not attempt here to perform a post-Newtonian expansion. This is a 
tedious task and not even well established for a theory which exhibits violation 
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of the EEP. Fortunately, this is not needed anyway for our purposes. Keeping 
just first order terms and applying the standard gauge conditions 

hi - \k, = , K,, - \h% = \hl, , (14) 
the time-time component of R^i,, for instance, reduces to 

Roo = -V=^(/ioo), (15) 

as is well known, where = dijdidj. Focus now on the rhs of eq. (12) and 
consider for simplicity just the second term in the square brackets. If EEP test 
are not to show sever violations, A should be chosen in such a way, so as 

2KXLm < 1, (16) 

since (1 + 2kAL^)~^ acts as a coefficient of T^i,. Assuming now, as usual, an 
almost static matter configuration, so that time derivative can be considered 

of higher order, one straightforwardly sees that the second term in the square 
brackets at the rhs of eq. (12) will lead to a term (amongst others) proportional 
to nW'^Ljn. I.e. 

V'(/ioo) oc fi:AV'L„ + E.S.T., (17) 
where E.S.T stand for extra source terms (not necessarily of higher order). 

The mere form of eq. (17) is enough to derive conclusions, and this is why 
a more detailed calculation is not needed: The time-time component of the 
weak field metric clearly has an algebraic contribution from L^. This implies 
that the metric is not plainly an integral over the sources, as in GR, but it 
also has a direct algebraic dependence on the matter field and their derivative 
(included in Lm). Obviously, whether or not Solar System test will be satisfied, 
or even whether the Newtonian limit itself could be recovered, would have to 
depend on the special characteristic of the matter, such as the energy density 
for instance. Additionally, any discontinuities or delta function in the matter 
could lead to unavoidable singularities! Note also, that tuning A in order to 
suppress this term cannot really ameliorate things here, since the problem 
does not really lie on the magnitude of this problematic term, but on its very 
presence (due to the possibility of exhibiting discontinuities etc.). 

Of course, as already mentioned, this is just one example of the problems 
that could be caused by the unusual differential structure of eq. (6) which we 
used for demonstrative purposes. It should be clear by now, even just based 
on the comparison with Palatini f{R) gravity and eq. (9), that the theory 
described by the action (5) cannot constitute a viable alternative to GR, for 
essentially the reasons described in Refs. [17,18,20] for Palatini f{R) gravity. 
As it has been already stated in [17] all of these issues do not constitute a 
special characteristic of the latter, nor are they specific (or generic) to the 
Palatini variation. They are just straightforward shortcomings of theories in 
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which gravity does not act in a cumulative way, i. e. theories in which the metric 
is not given as an integral but is also algebraically related to the sources and 
their derivatives, and therefore is influenced by local properties of the matter. 
Obviously, all theories in which the differential order of the matter field is 
higher of that of the metric in the field equations have this characteristic. 

Another way to understand this problem is to observe that in action (5) one 
has essentially added, with respect to GR, a strong coupling between the scalar 
curvature and matter, without however adding any extra dynamics. This is 
exactly the case for Palatini f{R) gravity as it becomes clear when the latter is 
seen as an ujq = —3/2 Brans-Dicke theory: The (auxiliary) Brans-Dickc scalar 
essentially introduces a strong coupling between gravity and matter (which, 
however, does not lead to violation of the equivalence principle) [8]. 

Notice that once any of the function /i 2 in action (2) are non linear in R, 
the field equations (3) immediately become fourth order in the metric. Now, 
even though the strong coupling still exists, more dynamics are introduced. 
The field equations are still third order PDEs in the matter fields, but now 
the differential order in the metric exceed the differential order in the matter 
fields, and therefore gravity becomes cumulative again. This implies that there 
is no reason to expect that the issues discussed here will persist. However, a 
more detailed study of the phenomenology of such theories is needed to arrive 
to definite conclusions. 

Finally, it is also worth pointing out that the action (5) constitutes another 
example case where the dynamics of theory arc not apparent in the action, 
due to the fact that the chosen representation of the theory is not helpful in 
this aspect (see [21] for a general discussion on theories and representations). 
Indeed, judging for action (5), one could easily be tricked to think that it will 
lead to a theory which is second order in the metric and only first order in 
the matter fields (as there arc only first derivatives of the latter). However, 
R includes second order derivatives of the metric. When R is not coupled to 
other fields (and has no self interactions), in a metric variation, these second 
derivative terms conveniently lead to a surface term. However, when R is 
present in coupling terms, by parts integration to arrive to total divergencies 
introduces extra derivatives of the fields R is coupled to; in this case this is 
the matter Lagrangian L^. 

To summarize, it was argued that a gravitational theory which, in addition 
to the Einstein-Hilbert term and the matter lagrangian includes also a direct 
linear coupling between R and the matter action cannot be viable. This has 
two important consequences: Firstly, if theories with coupling between R and 
the matter are to be explored further, cither the coupling or the gravitational 
action have to be taken to be non-linear (with all the complications this might 
bring). Secondly, a linear action with a linear coupling cannot be used even 
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as an approximation of a more complicated one for specific applications (as, 
for instance, in [16]), as it leads to spurious phenomenologjo and is in this 
sense unique among the general class of theories described by action (2) for 
arbitrary choices of /i^2- 
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